Self— consistent theory of the long— range order in solid solutions 
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On the basis of the assumption that atoms play a role of effective Fermions at lattice distribution, 
the study of the long-range ordering is shown to be reduced to self-consistent consideration of single 
and collective excitations being relevant to the space distribution of atoms and Fourier transform 
of such distribution, respectively. A diagram method advanced allows to elaborate complete ther- 
modynamic picture of the long-range ordering of the arbitrary compositional solid solution. The 
long-range order parameter is found for different chemical potentials of the components to obtain 
a scope of ordering solid solutions according to relation between degree of the chemical affinity of 
the components and mixing energy. The boundary composition of the ordering phase AB n is deter- 
mined as a function of the chemical potentials of the components and concentrations of impurities 
and defects. Temperature-compositional dependencies of the order parameter and the sublattice 
difference of the chemical potentials are determined explicitly. Polarization effects and passing out 
of the compositional domain 0.3f8 < C < 0.682 is shown make for transformation of the second 
order phase transition into the first one. The hydrodynamic behavior of the system is presented by a 
reactive mode being result of the interference of condensate and fluctuation components of collective 
excitations. The dispersion law of this mode is displayed experimentally as the Zener peak of the 
internal friction whose frequency and wave number decay monotonically with temperature increase 
and phase velocity has a maximum at intermediate temperatures in ordering domain. The polar- 
ization effects are shown to be relevant to the static component of Green function, the Goldstone 
mode of the symmetry restoration is represented by the instant vertex function. 

PACS numbers: 05.50.+q, 64.60.Cn, 61.66.Dk 



I. INTRODUCTION 

In spite of long history, the problem of long-range or- 
der in solid solutions keeps constant interest both of the 
academic community and for metallurgy applications 

- H . Main peculiarity of contemporary theories along 
this line consists in using two marginal approaches be- 
ing based on consideration either of the space distribu- 
tion of atoms of different kinds or the Fourier transform 
of such distribution named as concentration waves. In 
other words, long-range ordering problem is addressed 
to separate considerations of single or collective excita- 
tions distributed over thermodynamic states. However, 
it is quite clear that above excitations are strong cou- 
pled. Thus, the problem appears to study systematically 
both single and collective excitations in self-consistent 
manner. That is main purpose of the present paper that 
is based on the methods proposed initially in quantum 
statistics (?]]. 

The paper is organized in the following manner. Sec- 
tion |n| is devoted to statements of used formalism whose 
key point is that at the lattice sites distribution atoms 
of a given component is subjected to a prohibition rule 
to be considered as effective Fermions. This allows us to 
use the well-known formalism of the second quantization 
with generic Hamiltonian in a form inherent in the su- 
perconductivity theory. In Section III , a diagram method 



is built up for self consistent describing both single and 
collective excitations. This permits to get general equa- 
tions for the required set of the Green functions. As is 
shown in Section |y| devoted to consideration of the single 
excitations, their separate description is obtained in the 
simplest way within the framework of anomalous quasi- 
mean value. First, within this approach, we reproduce 
shortly main results of standard mean-field procedure 
and then we use the Green function method which allows 
one to describe not only each of the types of excitations 
but take into account their coupling as well. Practical 
advantage of this method is a possibility to study in stan- 
dard manner solid solutions with arbitrary magnitudes of 
chemical potential related to both different components 
and sublattices of ordering structure. As a result, we 
elaborate complete thermodynamic picture of the long- 
range ordering of arbitrary solid solution. Consideration 
of collective excitations in Section |v| allows us to repro- 
duce main details of this picture and to take into account 
polarization effects transforming the second order phase 
transition into the first one. The hydrodynamic behav- 
ior of the system is shown to be presented by a reactive 
mode being result of the interference of condensate and 
fluctuation components of collective excitations the last 
of which is of diffusive type. Concluding Section 
tains a discussion of results obtained. 
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II. EFFECTIVE HAMILTONIAN 

We consider a binary solid solution A-B with com- 
ponents A and B, which interact with potential vfjL = 
v ab (ri — r m ) if components a, b = A, B are placed in lat- 
tice sites l,m with coordinates r;,r m , correspondingly. 
A components distribution over solution sites is given 
by microscopic occupation numbers nf,nf — 0,1 being 
obeyed to conservation conditions 

nf+nf = l, l = l,2,...,N; 

£,n?=JV 01 a = A,B (1) 

where N is a total number of the solid solution atoms, 
Na, Nb are the same for components A, B provided N — 
Na + Nb- With accounting Eqs. (|l|), the configurational 
Hamiltonian of the problem 

W = ^EE<™">™-E^< ( 2 ) 

Im ab al 

where \i a is chemical potential of the components a, is 
reduced to usual form 

H = Ho + Hi + H2', (3) 
Ho=kEi m v?J?-»BN, (4) 

Hx = Ei(ei H)n h e t = £ m «f ~ «£f ), 

fi = fiA-^B, ni=nf, (5) 
n 2 = I J2im wi m nin m , w lm = vfj + vf* - 2v^.(6) 

Hereafter, the only A-component occupation number 
ni = nf will be used. 

The corner stone of our approach is that, in sense of 
the lattice sites distribution, atoms of a given compo- 
nent is subjected to a prohibition rule ni — or m = 1 
to be considered as effective Fermions distributed over 
"states" I ||. Thus, it is convenient to use the well- 
known formalism of the second quantization. With this 



aim, we represent the occupation numbers in the usual 
form ni = at at where creation and annihilation opera- 
tors af,ai are subjected to the anticommutation rules 
{a;,a^} = Si m (as usual, braces denote anticommuta- 
tor). Moreover, the ordering process causes lattice split- 
ting into two sublattices which will be labelled by Greek 
indexes or 1. As a result, Hamiltonian H = Ti — Ho 
counted out the non-essential constant Ho takes the stan- 
dard second quantization form: 

H = ^Ylei-fiala^aia + ^^wirna^a+arnoan. (7) 

la Irn 

Here, the anticommutation properties and natural con- 
dition wu = (a self-action is absent) take into ac- 
count, the expressions for a Fermion effective energy si, 
its chemical potential /i Q (being dependent on sublattice 
number a) and an effective interaction potential wi m are 
defined in Eqs. (H), (||). The key point is that atoms A 
on their own sublattice a = and on outsider one a = I , 
which chemical potentials are subjected to the condition 
fio < £i < /ii, are considered in similar manner as particle 
and hole in semiconductors. On the other hand, we are 
kept in last term the only contribution of the sublattices 
that causes the ordering process § . 

Sometimes, it is convenient to use instead of the above 
site representation the wave one: 

^#E^- ikri > <^Sj24e^', (8) 

j / 



ek ^ |pE e ' eikri ' Wk ^|]r^™ eik(r! ~ rm) (9) 

l Im 

where summation is fulfilled over N/2 sites of a sublat- 
tice. Then, the expression (R) takes the form 



H = g E l £ ~ M Q |aL a ka + E M k-k'akl a ik+q,O a k'O a -k'+q,l- (10) 
ka kk'q 

Here, we take into account the lattice translational invariance condition according to which the definition (||) derives 
to a constant Fermion energy si = e so that — £^ko- 



III. DEVELOPMENT OF THE 
SELF-CONSISTENT SCHEME 

For describing both single and collective excitations, 
it is convenient to proceed from different (site or wave) 
representations, so that we shall first make recourse to 
the diagram method making no use of the explicit form 
of the appropriate Hamiltonian. This will permit to get 



the general view of equations for the required set of the 
Green functions. 

Description of the present system is ensured by the use 
of Matsubara's Green function 

G?£(t) = -(Ta la (t)a+ f} (0)) (11) 

where at (t) , a/ Q (t) are the creation and annihilation op- 
erators within the Heisenberg representation, t is the 
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imaginary time, the angular brackets mean averaging 
over configurational states of the solid solution (see ||), 
the other symbols are standard 0. Making use of the 
above anticommutation relations, it is easy to convince 
that the normal (diagonal) and anomalous (off-diagonal) 
components possess the following properties: 



as follows: 



Gml{~ t) — @lm (*) = Glm(t), 
G ml(t) = G lm(. t ) = F lm{t)- 



(12) 



Thanks to this, it is sufficient to examine the behavior 
of only two components — the normal Gi m (t) and the 
anomalous Fi m (t) being shown in Figures la, b, respec- 
tively. Further, when making general computations, it is 
also convenient to use the matrix representation of the 
type Gim (t) where the cap signifies the exhaustive search 
for the a and (3 indices. 



FIG. 1: Diagram representation of single and collective exci- 
tations 

In addition to the matrix nature of the Green func- 
tion itself, the presence of the sublattices a, (3 = 0, 1 
results in matrix structure of the four-tail vertices of the 
effective interaction depicted in Figure lc (here, solid 
line corresponds to a solution component on own sub- 
lattice (a, (3 = 0) and on another's one (a,/3 = 1), the 
dashed line is respective for the mixture potential wi m ). 
Since the ordering is caused only by the coupling be- 
tween components belonging to different sublattices, for 
fixed indices !, rawe shall consider as non-zero only the 
components tuoi.oi = w ao , toio.io = w n> «>oi,io = w i, 
11)10,01 = wio producing the w matrix of the second rank 
(here, one means that interaction Wys,a@ transforms the 
sublattices a/3 into 7<5). The V matrix (see Figure Id) of 
the vertex function is structured in the similar manner. 
However, while the w matrix of the bare interaction is 
evidently diagonal, the complete vertex T has, as will be 
seen below, all nonzero-components. 

In the diagram representation the Dyson matrix equa- 
tion takes an ordinary form depicted in Figure le where 
the thick line corresponds to the exact Green function, 
and the thin line does to the bare one (the matrix of the 
latter is diagonal). The self-energy function (Figure If) 
is expressed by the equation illustrated in Figure lg. As 
a result, the problem reduces to the self-consistent de- 
termination of the vertex function T. It can be shown in 
standard manner [Q that, at the expansion in terms of w, 
the main contribution is made by the terms containing 
the polarizer II (see Figure lh) the matrix components 
Tl a p of which are determined in a similar way to the w a $ 
and T a p. Then the appropriate series reduces to the lad- 
der sequence which can be represented in the form of the 
Bethe-Salpeter equation shown in Figure li. It closes the 
system of equations for the self-consistent description of 
the ordering process in solid solution. 

In the analytic representation, this system is written 



G 1 (uj s ) = iui s 6 + Q2-S) 

S a /3(i) = Gf3 a (— t)r a p(t), 

f-\n s ) = w- 1 -n(fis), 
n Q/3 (t) = [G a0 {t)f. 



(13) 
(14) 
(15) 
(16) 



Here, the frequencies u s = ir(2s + 1)T, Qs = 2nST of 
the single and collective excitations, respectively, are in- 
herent in the Fermi and Bose particles, T is tempera- 
ture in energy units, s, S = 0, ±1, . . . are integers; 5 
is a unit matrix, e is a diagonal matrix with the ele- 
ments e 00 = — e 11 = £, £t is the same with the elements 
/j 00 = /io, /i 11 = — /Hi (e, jit a , a = 1, 2 are the bare energy 
and the sublattice chemical potentials of a Fermion) . The 
matrix structure of the Green function ( |l3| ) is stated on 
the Fermi conditions (|l2|). With derivation of the equa- 
tion (|l5|) for the vertex function T, it is generally taken 
that the bare potential reduces to the constant w. A 
distinctive feature of the system obtained consists in the 
fact that the explicit expressions ( |l3] ) , (|l5|) for the Green 
functions G, T of the single and collective excitations, re- 
spectively, are obtained in the frequency representation, 
while the expressions (|l4j) and ([l6]) for the self-energy 
function S and polarizer n require the application of the 
time representation. As far as the site and wave represen- 
tations are concerned, their choice depends on the type 
of excitations. 



IV. SINGLE EXCITATIONS 

If we are interested only in the behavior of the sin- 
gle excitations, their description is obtained in the sim- 
plest way within the framework of the quasi-mean value 
method ||. Applying the standard procedure, it can be 
shown on the basis of the expression ^ that in the limit 
N — * oo the behavior of the system is asymptotically 
defined by the approximating Hamiltonian which takes, 
in the self-consistent field approximation, the following 
form: 



H 



ef 



la I 



Anomalous quasi-mean values 

Vi = ( a to a u) ' Vi = (an«io) 



(17) 



(18) 



determining the amplitude of the effective Fermion trans- 
fer from one sublattice to another represent the (local) 
order parameter which value gives a gap 

= -- ^ WirnVm, A,* = - - ^ W lm T]* m . (19) 



Here, multipliers 1/2 in front of sums appear to take sum- 
mation over whole lattice sites (but not over sublattice 
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ones). In case of the long-range ordering, the order pa- 
rameter, given Eq. (|l8|), becomes constant (77; — > 77) and, 
within the usual mean-field approximation, Eqs. jl9| ) 
give the simple expression for the gap: 



A W 

= ~2 r1. 



W 



Wir. 



-zw 



(20) 



where the last equality is taken within approximation of 
interaction of near neighbors which mixing potential is w 
and number is z. For case of the short-range ordering, 



the correlation parameter ( jig ) slows down with I very 
fast. 



Diagonalization of Hamiltonian (17) is achieved by 
means of the transformations 



ai+ 
am - 



-- uia n 



v i a io , 



= uiaio 

u i a i+ 



man] 
■ viai- 



where, by virtue of the anticommutation rule, u 
It results in giving 



(21) 



= 1. 



J 



n ef = u 



6 t (a~l_ai- 



a T+ a i+) 



J_ai- -a l+ ai+) 



(22) 



where we are restricted ourselves with the case of a real 
order parameter (rji = 77^ ) and one denotes: 

u = fju- |Ei f i. V>= hiv-i - mo); 

ei = yftfi - ef + Af , /I = + m); (23) 



The quantity U represents the energy of the ground state 
whose wave function is 



|*o) =]J{ui +via+a m ) |0) 



(24) 



where |0) is the wave function of the Fermi vacuum. In 
general case of arbitrary concentration C of solid solu- 
tion, excitations defined by the operators ai+, ai- are 
non-coincident. Postponing this case below, let us con- 
sider at first the simplest case of stoichiometric solution 
(C = 1/2) where the coincidence of the behavior of ele- 
mentary excitations in the both states corresponding to 
sublattices a = 0, 1 is observed («;+ = = a). It is 
easy to see that in this case Jl = 0. 

The equation for the gap is obtained by substitution 
of the equalities (21) in the definitions (18), (|l9|): 



2 ^ 



Wlr. 



A, 

Err, 



tanh ■ 



2T 



A 2 . 

m 

(25) 
I 



Here, we have used Fermi distribution for excitations 
number 



vi = (at on) = [1 + exp(ei/T)]" 



(26) 



In the case of the long-range ordering, when dependence 
of the gap A; — > A on the site number I disappears, we 
obtain the usual equation within the mean-field approx- 
imation: 



2c 



tanh — = — , e = \/ (a ■ 
2T W v y 



W 



E 



Wlr, 



_ (27) 

Numerical solution of integral equation (E3) allows us to 
take into account effects of short-range correlations as 
well, for which the parameter A; is a strong function of 
the site number I. 

Concentrations C/ Q = (a^flfa) of A-type component 
on the sublattices a = 0, 1 is determined by Eqs. (|2l]), 

dD- 



it 



— tanh — — 
2T 



(28) 



On the other hand, the difference fi of the chemical po- 
tentials on the sublattices is given by the condition |J 



dn/djl= -NC\ 



(29) 



Making use of Eqs. (^2|) gives 

n = (u - n^ + tJ^H 1 - n) 



and thermodynamic equalities (p9|) arrive at trivial re- 



(30) 
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suit C — 1/2 for mean concentration. What about the 
entropy S = —dfl/dT, thermodynamic potential (|3C| ) de- 
rives to usual configuration form 



S = -J2 + C 1 - vt)Hl - 



(31) 



With transition to the ordered state, the jump of the 
thermodynamic potential AF = F orc i — Fdi SO rd is deter- 
mined by the formula 



AF 



2W ^ 



(32) 



following from the equalities dF/dAi = (dHef /dAi) and 
( p2| ) , (|2^) , ( p6| ) , (|7|) . According to ( |32] ) during ordering, 
as it usually is during transition to the low-symmetry 
phase, value F decreases. 

The elementary excitation energy £ ex = ^2 e l u l IS at 
C = 1/2 of the form: 



E 



V(m-£) 2 + A 



; 1 + exp 



^(t-ey + Af/T 



(33) 



Previously, we applied the mean-field method bringing 



in shortest ways to the description of single excitations. 
Naturally, the illustrated results can also be obtained by 
means of the Green function method which allows for 
determining the Fermi function G(lo s ) from the system 
( |l3| ) — (^6|)- Keeping in mind the fact that the Green 
procedure allows for describing not only each of the types 
of excitations but for taking into account their coupling, 
we shall demonstrate at first the way it reproduces the 
results of the mean-field method. 

In the context of the above used site representation, the 
functions G, e and E in the equation ( |l3| ) and functions 
r, w and II in the equation ([l5]) should be corresponded 
the lattice index Z, and factorization of the appropriate 
contribution ( |l4| ) and ( |l6| ) is achieved in transition to the 
wave representation. 

To define the explicit form of the Green function (|l^) , 
we shall consider that, as in the case of the theory of 
superconductivity, the self-energy function E/ acquires, 
under ordering, the off-diagonal components Eoi ; = 
E^q ; = A/ corresponding to the gap A; in the energy 

spectrum of the single excitations 0. The matrix Gi 
inversion leads in the case C = 1/2 to the following ex- 
pressions: 



J 



G l (tu s ) = - 



\uj s + ((j, - e) 



Fi(lo s ) = — 



= J(e-/2) 2 +Af, uj s =7r(l + 2s)T, s = 0,±L... (34) 



r 



With accounting the sum rules 



T E 



s— — oo 



tanh — 
2e IT 



. E 



(35) 



being applicable for some real e, substitution of Eqs. (|34|) 
into Eq. ([l4]) for the off-diagonal components, where the 
vertex T/ is reduced to convolution operator of the bare 
potential wi = J2 m Wlm > arrives at the self-consistency 
equation which, as would be expected, coincides with Eq. 
( p5j ) for the gap A; . It is key point that the diagrammatic 
approach allows to approve the self-consistency solution 
of Eq. (25) by means of using some approximation for 



the vertex T/ (see Section [v]) . Dependence of the effec- 
tive chemical potential Jl on the difference of the lattice 
concentrations C/ Q follows from the relation 



C, 



; {C 11 -C iq ) = -TY j Gi{u 



(36) 



which, with the application of the first one from the 
equalities (53), leads to result 



Q = ^£tanhiL 
2e; 2T 
following from Eq. (Eq) as well. 



Polarizer (nff) has the following Fourier-form: 



(37) 



oo 1 



N 

S— — OQ m 



(38) 

Substituting the Green functions (p4) here, upon sum- 
ming over s, we arrive at the expressions 
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n 00 (0) = n 11 (0) = -27v^er 1 tanh^+AA 2 , n 01 (0) = n 10 (0) = -AA 2 ; (39) 

I 

A = jn J2 e i 3 (tanh — ^cosh 2 -§f) , 
i 



where, in view of the macroscopic equivalence of sites, 
there is no dependence on their number I. Making use 
of Eq. ([38]), we find the inverse vertex function ( |l5|) re- 
sponsible for the behavior of the collective excitations 
(see Section |y|) . 

The system of equations (Q), (pD, (|25|), (|29|) and 



that offers the self-consistent description of the single ex- 
citations affords finding the gap A, the order parameter 
r\ and the sublattice difference of concentrations C by 
the assigned values of thermodynamic parameters (the 
temperature T and the chemical potential p), as well as 
microscopic parameters (the value of the mixing poten- 
tial wi m and effective energy e given by Eqs. (||), (||)). 
Numerical solution of these equations for the long-range 
order parameter results in the temperature dependencies 
displayed in Figure 2. From this figure we notice that 



FIG. 2: Temperature dependencies of the long-range or- 
der parameter at different magnitudes of chemical poten- 
tials (curves top-down correspond to values [fi — s)/Aq = 
0.0, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9. Concentration 
C= 0.5) 

the parameter 7/ monotonically decays as the tempera- 
ture T and the difference \p — e\ increase. According to 
Eqs. (|27|), (|37|), at long-range ordering the sublattice 
concentration difference is 



to be determined by the mean chemical potential but 
being non-dependent on temperature. Usually, one sup- 
poses p = e and, as consequence, C = (i.e., Co = C\). 



C 



I* 



(40) 



The ground system state defined by the function ( pi ) 
is achieved at T = where site dependence is miss- 
ing and excitations are absent [y = 0). At p = e the 
gap width and order parameter take maximum values 
Aq = i;W = —\ J2m w im, f] = 1, variation of the ground 
state energy, being equal to the decreasing thermody- 
namic potential, takes maximum value |A.F| = -^A ; 
respectively, the sublattice concentrations Co, C\ co- 
incide. With finite d ifference p — e ^ 0, t here are 
A = VAg-Oi-e) 2 , V = ^1 - (/2 - e)VA§, AF = 
-f A [1 - (M - e) 2 /Al] ,C = (fi- e)/2A . Thus, as 
the sublattice chemical potentials fiQ = \i\ = /i increases 
from value e to e + A , the gap width, order parame- 
ter and absolute value of the system free energy change 
monotonically decay going into zero at the boundary 
value p c = e + Aq. The sublattice concentrations dif- 
ference varies here linearly from to 1/2. 

As the temperature arises, originated are the el- 
ementary excitations the number of which v w 
exp(— Ao/T) <C 1. With a precision of the first power 
terms over v <C 1, we obtain: 



A = A <i/l- 



An 



1 - 2 



exp (-A 



<0 



AF 



fA ?7 2 , £ cx -7VA exp(-A2 



(41) 



As it is known, the exponential character is inherent in the low-temperature dependencies usually observed under 
phase transitions. 

Near critical temperature T c = W/A the ordering process can be presented in analytical form within the mean-field 
approximation that is stated on making use of Eq. (p7|). Within domain < (T c — T)/T c -C 1, we obtain: 



A = 2T C VM, ry = V%i, t 



T a —T 



\w, w = -J2 m w lm , 



AF = -%NT c t, £ cx = NT c V3i 



(42) 



where we put ji = e. 
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Now, we consider general case of non-stoichiometric solid solution which the (mean) concentration C ^ 1/2 and 
the difference fx of the sublattice chemical potentials is determined by equations (cf. Eq. (^)) 



dn/djl= -NC; 



fi = -f7i- i£ Q + T£m(l-^+)+T£ m (l-^-), (43) 



i 



vt± = {l + exp(^g2)} 

following from Eqs. (f29|), ( |22] ) and (p3|). Respectively, inversion of matrix ( |l3| ) gives instead of Eqs. ( |34| ) 

G°°K) = [iw, - ( m - e)] /V h G} x {uj s ) = [kj, + (no - e)} /V u 
Gf{u s ) = Ai/V h Gj°(co s ) = Af/Vf, V t {uj s ) = [iuj s - (e, + /*)] K + ( CJ - £)] , (44) 

e i = v / (M-^) 2 + l A d 2 I M=|(i"o+Mi)) M= |(a*i -Mo)) 
cj s = tt(1 + 2s)T, s = 0,±l,... 



With accounting Eqs. ( |3 5|) and Eq. (14) with T substituted by w, we obtain the self-consistency condition 

-7 E — ( tanh + tanh e -^) = A < C 45 ) 



in x ' 

replacing Eq. (p5|). Respectively, from equalities C a = — T^2 s G aa (LJ s ) one follows: 

= e -(jj- e) e + Ti _ e +(M- g ) 
u 4e 2T 4e 2T ' 

= e+(/2- g ) e + M _ e-jfl-e) e-£ 
4e 2T 4e 2T 



(46) 



In case of the long-range ordering, when the parameters A; — > A, e; — > e do not depend on the site number Z, the 
relations (|45|), (46) arrive at Eq. (pp|), as follows. On the other hand, thermodynamic equality ( f43| ) derive to equation 

22-1= -T N+-^-) = — T (tanh tanh (47) 

AT ^ V + ' 2N ^ \ 2T 2T J V ' 

i i v 7 

where the relation 1 — 2vi± = tanh £ ' ? ^ is taken into account. In difference of Eq. (p9|), this equation is non-trivial 
and accompanied with Eqs. (f45|), (^6|) determines quantities A;, Co, Ci, /I at given thermodynamic parameters T, 
C, /I and microscopic ones wi m and e. In similar manner, the excitation energy is determined by the expressions (cf. 
Eq. @) 



P \^ / , \ a f-, , cosh(e ; /2r) ) 1 
£ cx = ^ e; (, i+ + .*_) = A ^ |1 + CQsh(/7/2r) ) 



(48) 



at fjb = e. The entropy 5 = —dVL/dT is reduces to the configuration form type of Eq. (|3l|): 

S = - ih+ In n+ + (l - n+) hi(l - v l+ )] + [n- In n- + (l - ^-) Ml - n-)]} • (49) 
i 

In the simplest case p = e equations (]45|) , ( |47| ) govern the long-range ordering to be reduced to the form 

ri + fh n — fh , 

tanh-^ htanh — = 2ry, (50) 

tanh _ tanh = 2(2C-1) (51) 

26 26 v ; v ; 



I 

where dimensionless values of the sublattice difference of the chemical potentials fh = /js/Aq and temperature 6 
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T/T c are introduced. Analytical consideration of this 
system is possible in the limits 9< 1 and t|<l only. 
In the first case, with accuracy of the first non-vanished 
terms, equations (|50|), ( fill) arrive at the expressions 

^ = r/0 -exp(-^o), (52) 
Jl = A i] Q + aT (53) 
where parameters 770, a are determined by relations 



1 , a 

Vo = 2 ( * 2 



(l-Cf 



2(C- 



3)2- 



(54) 



The phase transition curve where 

, TO 1 1 + S 2 

rj = 0, cosh- 



e ve 1 

is determined by relations 



5 2 '' 



5 = 2C-1 (55) 



2C* - 1 



1%. 



T 



T = T r 



1-S 2 
1 + S 2 



In the vicinity of the transition curve T C (S) one has 
r, 2 = A(S)t(S), t(6) = ^^; 



cosn 2T — 1(52 



A(8)8 2 



-l=jr*(*) 

where factor j4(5) is determined by the equality 



A~ 



2S 2 



1 + S 2 



1 



1 + 5 2 



1 



l-<5 2 



(56) 

(57) 
(58) 

(59) 



With deviation of the stoichiometric composition this fac- 
tor decays as A ~ \/3(l — 2S 2 ), S 2 <C 1 taking zeroth mag- 
nitude at <5o = 0.364 (Co = 0.682). Physically, this means 
that with passing out of the domain 0.318 < C < 0.682 
the second order phase transition is transformed into the 
first one being close to the second. 



FIG. 3: Temperature dependencies of the long-range order 
parameter (a) and the sublattice difference of the chemical 
potentials (b) at different concentrations C being pointed 
out near corresponding curves. (The dashed line bounds the 
abruption region) 

Explicit forms of the dependencies 77 (T, C), J1(T, C) in 
the case p, — e = are depicted in Figures 3, 4. The 
order parameter 77 decays monotonically with increase 
of both temperature and concentration (see Figures 3a, 
4a) , whereas behavior of the sublattice difference of the 
chemical potentials Jl is much more complicated: when 
the temperature grows, it decreases near the stoichio- 
metric composition C — 0.5 passing to increasing regime 
with approaching to the composition C = 1 (Figure 3b); 
analogous behavior is observed with deviation of the sto- 
ichiometric composition at low and pre-critical tempera- 
tures (Figure 4b) . Principle important is the interruption 



FIG. 4: Compositional dependencies of the long-range order 
parameter (a) and the sublattice difference of the chemical 
potentials (b) at different temperatures (the magnitudes of 
the latter related to critical value are pointed out near cor- 
responding curves). The dashed line bounds the abruption 
region 



of the sublattice difference of the chemical potentials Jl 
at the curve of the phase transition that takes maximal 
values within intermediate ordering region. 

If the system is subjected to external influence, the 
parameter to = (/2 — s)/Aq becomes non-zeroth and or- 
dering process is described by the system (dO), (51 ) where 



the order parameter 77 should be replaced by the renor- 
malized one \Jr\ 2 + to 2 . It is easy to convince oneself the 
growth of the influence parameter to causes suppress- 
ing ordering process to shrink the long-range domain as 
shown in Figure 5a. At that, the boundary composition 
being the solubility limit of the ordering phase is deter- 
mined by the parameter of external influence in following 
manner: 



C = l 



y/m(l — to) — TO 
2(1 - 2m) ' 



(60) 



According to related plot in Figure 5b the solubility con- 
centration decays monotonically with arising external in- 
fluence to vary anomalously fast near boundary magni- 
tudes C = 1 and C = 0.5. 

FIG. 5: a — The phase diagram determining long-range order 
(LRO) and disorder (DO) domains at different magnitudes 
of chemical potentials (curves top-down correspond to values 
{p-e)/A = 0.0, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9). b- 
The solubility limit of the ordering phase versus the parameter 
of external influence 



COLLECTIVE EXCITATIONS 



As is seen from Sections [II| and [V , the ordering process 
results from the fact that, as the temperature decreases, 
the effective Fermions form coupled pairs corresponding 
to the collective excitations of the Bose type. If the num- 
ber of such pairs Af = \Nr) makes up the finite part 
with respect to the total number of sites N, the ordering 
process assumes the macroscopic character and is deter- 
mined by the long-range order parameter 7/ = A/Aq. 
So, in the representation of single excitations, the solid 
solution ordering is reflected through the appearance of 
the off-diagonal components A in the matrix of the self- 
energy function E. 

Let us show now how the ordering manifests itself in 
the representation of the collective excitations. Appli- 
cation of the mean-field method, therewith, turns to be 
insufficient and a recourse should be made to the self 



consistent scheme developed in Section III to take into 
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account the behavior of the Green functions of both single 
and collective excitations. To study the latter we start 
with consideration of the Fourier representation T(K, fi) 
of the vertex function acquiring under ordering a conden- 
sate component r\)(K, SI) accompanied by a fluctuation 



function T'(K, O). Their principle difference is in hydro- 
dynamic behavior in the limit K, 57 — > where the latter 
tends to zero, whereas the former takes a finite magnitude 
To 7^ 0. For determination of the latter we substitute the 
polarizer ( ]38| ) into (|T5|): 



J 



f-i 
1 o 



B-A\A\ 2 AA 2 
A(A*) 2 B~A\A\ 2 



B 



-W~ 



1 . e 
— tanh — — 
2e 2T 



(61) 



r 



where account is taken of the diagonal structure of the in- 
teraction matrix w, the parameter A is determined by the 
last equality ( |39| ) (in this Section, the gap A is supposed 
to be complex). It is principally important the parame- 
ter B — > with ordering, in accordance with condition 
Eq. (H) for the gap |A|. 

Inversion of the matrix ( |6l| ) arrives at the off-diagonal 
component of the condensate vertex: 



1 o 



-A 



itanh^ - J-cosh" 2 ^ 
(l-Itanh§)(l->sh- 2 §) 



(62) 



Near the critical point (T — T c <^ T c ) where 

,2 



-.01 



2 -A V- 

3 °{1~T/T C ) 2 ' 



(63) 



we can neglect a frequency dispersion to put £° (fi) ~ 
£ 01 (0), r 01 (q) ~ r 01 (0) = T^ 1 in Fourier transform of 



(64) 



the equality (|14|) that takes the form 



1 . T] 

— tanh — 
r\ 6 



,Ao 

1 



As a result, the collective excitations are characterized 
by the relations 



V3 1 



T 



r: 



01 



-2A ; T-T C «T C (65) 



the first of which is principally different of the corner pe- 
culiarity ( |57j ) inherent in the single excitations. This dif- 
ference is expressed physically in that the collective exci- 
tations polarize ordering system to transform the second 
order phase transition into the first one. 

Now, we are in position to study the collective excita- 
tions themselves. Their Green function is determined by 
the fluctuation component r'(K, f2) of the vertex func- 
tion according to the relation 



^(K,fi) = -— r'(K,0). 

27T1 



(66) 



The simplest way to find this function in hydrodynamic 
limit K, fl — > is to apply the method developed in 
Refs. pT| , |12|] . Its main point consists in the fact 
that instead of the equation G -1 = Gq 1 — Ao<5, for the 

Green function G, use is made of its analog (g^ 

£-1 



- Aq5 for the averaged function ^Gy and locator 
Go)- In addition, inserted is the effective interactor 



U = AqS + Aq ^Gy which contains the term in the first 
order over Fermion energy. The Dyson equation replac- 
ing (|l3l ) therewith assumes the form U^ 1 — A^S — a. 
The two-particle Green function 



<p(E; K, 0) = ]T (G R (k+, kV; E + fi)G A (k_, k'_; E)) 



k,k' 



appears, in the ladder approximation, as (compare with ( pL 5| 

1 



<p(E;K,Q) = - 



2m 



7 



1 -^C/ R (k + ;£; + f7)C/ A (k_;i?) 



(67) 



(68) 



where k± — k ± K/2, 7 is an irreducible four-tail vertex, indices R and A of the retarded and advanced functions 
correspond to the selection of the components a = 0, 1 in various sublattices. Hence, taking into account the Dyson 
equation and Ward identity [Q , we have 



< (E + fl) - at (E) = 7 £ [C/ R (k' + ; E + fi) - C/ A (k'_ ; E)] - 7 fi 



(69) 
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where the irreducible vertex 7 has, contrary to the 7, two 
coinciding tails appropriate to the same sites. Then, we 
come to the conventional expression for the fluctuation 
component of the two-particle Green function 

^ K ^—nTimJW (70) 

where x(K) is susceptibility taking in hydrodynamic 
limit K — * the thermodynamic value %, D(K, SI) is the 
dispersing diffusion coefficient. As a result of the above, 
for the fluctuation component of the vertex function we 



obtain || 

(f< K .«>r^(- in r K2 ia + ° CK >> 

where the matrix structure reflects the presence of two 
poles SI = Ti-DK 2 . 

The sum of the expression ( |6l| ) with B = and Eq. 
( [fl| ) produces, after the matrix inversion, the complete 
Green function of collective excitations: 



0( KSl )--!S^(-to- i (hSKo-DK*) -i(A/\A\) 2 SK \ f72l 
<P^,U)- v(K,n) y -±(A*/\A\) 2 SK Si - i (±S^o - DK 2 ) J ' 

2?(K, fi) - f! 2 + ((1/2)5^0 - DK 2 f - (\/A)S 2 K 2 , (73) 

S 2 = Att X A\A\ 2 D, K 2 = 4t7 X A\A\ 2 /D. (74) 



Condition 2?(K, f2) = brings to the dispersion law 

SI = TDK^K 2 -K 2 (75) 

whose characteristic feature consists in a bell form with 
a slanting long-wave side (see Figure 6). Physically, the 
dependence (JT^) relates to the collective mode that is of a 
reactive nature in the long-wave region K < Kq limited 
by the wave number given Eq. ( f74|) and has a relaxation 
nature in the short-wave region K > Kq. At K ^> Kq, it 
transforms into the ordinary diffusion mode SI = =pi_DK 2 . 



FIG. 6: Law of collective mode dispersion (the solid line cor- 
responds to the real value of frequency Q, the dashed — imag- 
inary value, the dotted — diffusion mode) 

In this way, the self-consistent consideration of single 
and collective excitations leads to the conclusion that 
apart from the diffusion regime realized in the meso- 
scopic region K > Kq, possible in the system is quasi- 
oscillations characterized by the phase velocity S deter- 
mined by the first of equalities ([74]) . Obviously, the dis- 
persion law ( f75| ) acquiring the acoustic form in the limit- 
ing long-wave region K <g; Kq is observed as Zener peak 
in experiments with internal friction [ |j"3"| . This mode is 
characterized by the frequency 

Sl = SK - ^- ( - tanh 1 - icosh" 2 1) (76) 

whose temperature dependence at concentration C = 0.5 
is shown in Figure 7 with solid line. This frequency is seen 
to decay monotonically from the magnitude SIq — 7rx/Ag 



at T = to £1 = at T = T c varying near the critical 
point quadratically: 

Q ~ 2tt^- (l - |^ , T-T C «T C . (77) 

Here, we take into account polarization effects trans- 
forming the corner peculiarity ( p7j ) into the linear 
relation (|65|). Since the diffusion coefficient D = 



FIG. 7: Temperature dependencies of the characteristic 
frequency flo (solid line), phase velocity S (broken line) 
and boundary value of the wave number Ko (dotted line). 
Magnitudes Qo, S and Ko are measured in units 7rx/Ao, 
\j 7rx-Do/Ao and \J 7rx/Ao-Do, respectively. The height of dif- 
fusion barrier is put to be Q = 3T C . The insertion is built with 
accounting polarization effects given by Eqs. ( frl| ) — ( |79| ) 

Dq exp{— (Q/T)} is determined by the energy barrier Q 
to behave itself in a non-critical manner, the parabolic 
dependence (^) means the linear behavior for both the 
sound velocity and the boundary wave number: 

^V^S- 1 (£)('-£) (79) 

the former of which is exponentially smaller than the 
latter. With the temperature falling down in the domain 
T <C T c the sound velocity decreases and the boundary 
value of the wave number, on the contrary, infinitely in- 
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creases: 



K 



A A) 



■ex Pl -- 



Q_ 

2T 



■ exp 



(80) 
(81) 



A concentration deviation off the stoichiometric composi- 
tion causes order suppressing to decrease the values fio, 
Kq and S. It follows therefrom that the most prefer- 
able (in terms of detection of the zero-sound mode of 
collective excitations) is the temperature region that is 
situated just below the critical temperature T c . Explicit 
form of the temperature dependencies Qo(T), S(T) and 
Kq(T) are depicted in Figure 7 where main panel relates 
to the phase transition of the second order, whereas inser- 
tion takes into account polarization effects transforming 
it into the first order. 

According to Eqs. (Q — ( |rij ) the diagonal compo- 
nents of complete vertex function have the following hy- 
drodynamic form: 



r Q «(K,r! s ) = 27r X 



n 2 s + n DK 2 ' 



(82) 



n s = 27rST, S = 0,±1, 



Then, taking into account the sum rule jTo) (cf. Eqs. 



S=-oo 



Q 2 S + e 2 2e 



1 , e 
— coth — , 

2T' 



(83) 



for instant vertex function T aa (K) = TJ2 S T aa (K, fi s ) 
we obtain 



r-(K ) = ^/^coth^ g? . (84) 
y J 2 A V DK 2 2T v ; 



In the hydrodynamic regime K <C T / y/CloD this expres- 
sion arrives at the known Bogolyubov singularity H] 



r-(K)= 7 r^ 



(85) 



that keeps this form for off-diagonal components of the 
vertex function as well. 



VI. CONCLUSIONS 

Above consideration shows that within the single exci- 
tation representation the study of the solid solution or- 
dering can be achieved by analogy with the microscopic 
theory of superconductivity Q , Q . Description of the 
collective excitations requires their self-consistent con- 
sideration along with single ones. The hydrodynamic be- 
havior of the system is presented as the result of interfer- 
ence of the condensate and fluctuation components the 
last of which is of diffusive type. This interference results 



in appearance of the reactive mode corresponding to the 
zero-sound. 

It is naturally, the scheme developed allows one to re- 
produce complete picture of the solid solution ordering 
within mean-field approximation - - || . Setting aside 
well-known information along this approach, let us focus 
on main results of our consideration that are not achieved 
within mean-field scheme. 

First, we have found the temperature dependence of 
the long-range order parameter for different differences 
of the chemical potentials of the components fi = fiA^^B 
accounted from Fermion energy e (see Figure 2). Taking 
the latter as the origin, we obtain a scope of ordering 
solid solutions is determined by the condition 



~ Ms < A = W/2 



(86) 



fixing degree of the chemical affinity of the components 
with respect to the characteristic value W = — 7" ' w; m 
of mixing energy (^). In accordance with Eq. (|40|), it 
means the difference of sublattice concentrations 2C = 
C\ — Co could not take magnitudes more than one. 

The second of practical results allows us to determine 
the boundary composition of the ordering phase AB n : 



C = 



\J m(l — m) — m 
(n+ 1)(1 - 2m) : 



I^A - 

A ' 



(87) 



As distinct from equality ( p0| ) related to the simplest 
phase AB, here we consider the concentration domain 
< C < (n + l)- 1 instead of 0.5 < C < 1. Accord- 
ing to relation (pTfi, in perfectly ordering solution where 
fiA = fJ-B , the solubility limit is C — and grows anoma- 
lously fast with decreasing degree of the chemical affinity 
of the components with respect to the mixing energy. 
Moreover, if we insert impurities or defects with concen- 
tration x <^ 1, the parameter of external influence varies 
as 

m = mo + /3x (88) 

where the value mo is inherent in the proper solution 
AB n , (3 is constant taking both positive and negative 
signs. Thus, alloy doping arrives at lowering boundary 
concentration in the case (3 < and its growth otherwise. 

Principle advantage of our approach is a possibility to 
study solid solutions with arbitrary composition to elab- 
orate complete thermodynamic picture of the long-range 
ordering process. According to Figures 3, 4 the order pa- 
rameter decays monotonically with increase of both tem- 
perature and concentration, whereas the sublattice dif- 
ference of the chemical potentials suffers abruption with 
maximal value Jl ~ Ao/2 at phase transition related to 
intermediate compositions. This allows us to estimate 
the surface tension coefficient a ~ ji£,/v as follows: 



0.1 



(89) 



where £, v are characteristic values of the correlation 
length and the atom volume, respectively. Another phys- 
ical result is that the second order phase transition is 
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transformed into the first one with passing out of the 
compositional domain 0.318 < C < 0.682. However, a 
difference between these transitions is so weak to be ob- 
served in Figure 3a. 

Much more important is the influence of collective ex- 
citations whose polarization effects transform the phase 
transition to the first order always (see the first of Eqs. 
(|65|)). Related collective mode is of a reactive nature in 



the long-wave region limited by the wave number (74). 
The dispersion law (ffq) of this mode can be observed 
as the Zener peak of the internal friction. According to 
Figure 7 proper frequency and boundary wave number 
of this peak decay monotonically with temperature in- 
crease, whereas characteristic velocity has a maximum at 
intermediate temperatures in ordering domain. A devia- 
tion off the stoichiometric composition causes suppress- 



ing of the above zero-sound mode to be pronounced in 
the temperature region just below the critical tempera- 
ture. It is worthwhile to notice the polarization effects are 
relevant to the static (condensate) component of Green 
function ( f72| ) related to frequency = 0, whereas the 
Goldstone mode of the symmetry restoration is repre- 
sented by the instant vertex function (|85|) related to time 
t = 0. 
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